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ABSTRACT 



We present a fast method for reconstructing Differential Emission 
Measures (DEMs) using solar coronal data. On average, the method 
computes over 1000 DEMs per second for a sample active region ob- 
served by the Atmospheric Imaging Assembly (AIA) on the Solar 
Dynamics Observatory (SDO), and achieves reduced chi-squared of 
order unity with no negative emission in all but a few test cases. The 
high performance of this method is especially relevant in the context 
of AIA, which images of order one million solar pixels per second. 
This paper describes the method, analyzes its fidelity, compares its 
performance and results with other DEM methods, and applies it 
to an active region and loop observed by AIA and by the Extreme- 
ultraviolet Imaging Spectrometer (EIS) on Hinode. 

Keywords: Sun: corona - Methods: numerical - Techniques: image 
processing 



The Sun's outer atmosphere, the corona, is a hot (~ 10^ Kelvin), highly 
magnetized, dynamic plasma, with spatial scales ranging from the electron 
gyroradius of less than 1 cm to the solar radius of over 10^^ cm. As such, it is 
one of the most fruitful objects for the study of astrophysical plasmas. One of 
the most fascinating puzzles of the Corona is how it is heated to temperatures 
~ 1000 times those of the underlying photosphere. It is believed that this 
process is magnetic in nature, but the details (e.g., the spatial distribution 
of heating and whether it is continuous or impulsive) are not known. One 
important tool for understanding these details is the reconstruction of solar 
temperature distributions called 'Differential Emission Measures' (DEMs) from 
spectral images of the Sun. 



The DEM, £(T) characterizes coronal emission at a given temperaturq^l 
and is defined in terms of the densities, n(/), along the line-of-sight so that 



For most spectral lines; lines corresponding to forbidden transitions have 



1. Introduction 
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/ £(T)dT = J n^(l)dl (column emission convention). Observed intensities, /„, 
are given by integrating the instrument response, Rn{T), against the DEM being 
observed: 

4 = / Rn{T)8{T)dT (1) 



Algorithms which reconstruct these temperature distributions must be 
very fast if they are to process the volume of data produced by modern solar 
observatories, or resolve the dynamics and fine spatial scales believed to 
be involved in coronal heating . For example, matching the real-time data 



rate of AIA (jLemen et al.ll2012h , requires computation of over 10^ DEMs per 
second. There are a number of existing algorithms for reconstructing DEMs 
from solar image data, but they are far too slow to meet this requirern ent. 



One widely used method is the PINTofALE code (|Kashyap & Drake 



19981), which employs a Markov Chain Monte Carlo (MCMC) search 
taking seconds or minutes to compute a DEM. Another method, which 
is considered rela t ively fast and was recently applied to AIA data by 



Haimah & Kontail {iQl'^ . computes ~ 4 DEMs per second (an updated version 



with substantially unproved performance has recently been made available 
at http : / / www . astro . g la.ac.uk/~iain/ demreg/map/; we use this 
version in our comparisons). 

We present a fast, iterative, regularized method for inferring DEMs using 
data from solar imagers such AIA and EIS. With one thread on a 3.2GHz 
processor, it is able to compute well over 1000 DEMs per second for an example 
solar active region observed by AIA, and over 100 DEMs per second in our 
most difficult test cases. Moreover, we anticipate that with straightforward 
optimizations (e.g., conversion of computationally intensive portions of the 
code to C and parallelization), the performance of the code will be increased 
to ~ 10^ DEMs per second on a single workstation, sufficient to match AIA's 
real-time observing rate. 

This paper describes the method, analyzes its fidelity, compares its 
performance and results with other DEM methods, and applies it to example 



more complicated density sens itivity, but that is beyond t he scope of this paper. 
Interested readers may consult Idel Zanna & MasonI (|2005h . 
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solar data. We also touch on the limitations of the solar data to constrain the 
DEM being observed. It must be noted that the ability to recover the details of 
the original solar input DEM is limited, both because of the number of channels 
are limited and because the temperature response functions are broad, due 
to the width of the temperature-dependent emissivities of the spectral lines 
from which they are c onstructed. This i s discussed in detail, for instance, in 
Craig & Brownl (|1976l) : lTudge et al. 



2. DEM Algorithm 
2.1. First Pass 

We wish to infer the DEM, a continuous function, from a small number of 
observed intensities which are the result of the convolution of the DEM with 
the instrument response functions. This problem is inherently ill-posed due to 
the limited number of response functions and the loss of information resulting 
from the convolution (i.e., equation [T]) . To resolve the ambiguity, we impose 
additional constraints on the DEM solution. The first is that the DEM can be 
expressed as a linear combination of some set of basis elements (e.g., a set of 
narrow temperature bins), Bj, with coefficients ej (i.e., S{T) = ^jBj)- The 
integral equation ^ then becomes a matrix equation: 

/„ = j Rn{T)£{T)dT = Y,[ekj Rn{T)Bu{T)dT\ = J^^klnAnk, (2) 

k k 

Where 7„ are a set of normalization constants for the response functions, which 
we choose to be their squared integral, square root. 

If the number of basis elements is greater than number of instrument 
channels, the underco nstrained invers ion can be resolved by a Singular Value 



Decomposition (SVD; iPress et al.ll2002l) , which picks the coefficient vector with 



the smallest magnitude (i.e., least squared emission measure). This is a sort of 
smoothness constraint, since a solution which has emission concentrated into 
narrow peaks will have greater total squared EM than a solution which spreads 
the emission as broadly as the data will allow. Such a solution will also tend to 
reduce non-physical negative emission, since negative EM is likely to require 
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excess positive emission measure elsewhere in order to produce the observed 
intensities, resulting in relatively high total squared EM. 

If the number of basis elements is equal to the number of instrument 
channels, it is straightforward to invert equation |2] to find the DEM coefficients, 
Cfc. In that case, however, the quality of the inversion is highly sensitive to 
choice of basis. In particular, many choices of basis lead to an inversion matrix 
(i.e., Ajj^) that has large negative entries, causing large negative values in the 
resulting DEMs. 

The tendency of the inverse matrix to have large negative values is 
reduced, however, if Ajk is diagonally dominant (i.e., each diagonal entry of 
the matrix is the largest one in it's respective row/ column). In particular, if we 
choose the basis functions to be the instrument response functions themselves, 
square-normalized, we obtain a symmetric Aj^ matrix whose diagonal entries 
are unity, with all off-diagonal entries less than one. Remarkably, we find 
that this basis gives results identical to using a large number of narrow basis 
elements and an svd. It therefore satisfies the SVD's minimum squared 
emission measure constraint, and we use it to compute a first pass DEM as a 
starting points for our DEM inversion. It reduces the number of operations 
required for each pixel to only ~ A^^ operations with instrument channels, 
increasing the speed of the inversion. 

Even when using the response functions as a basis, however, we continue 
to compute A^j^ using an SVD. This allows us to enforce a minimum condition 
(~ 10^^, typicall}]^ on the inversion, ensuring that round-off error (from Aij 
being nearly singular) does not cause ringing. In effect, this combines highly 
similar channels (or linear combinations of channels, to be more precise) into a 
single channel rather than trying to fit them individually. 



^Appendix lAl demonstrates why this is the case 

■^This means the smallest singular value used to form the inverse matrix must 

1-12 



be at least 10 times the largest singular value. 
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2.2. Regularization 

Another constraint which DEM solutions must satisfy is a physically 
required positivity constraint. However, the first pass DEM inversion described 
in Section IZT] can produce negative coefficients and therefore negative EM. This 
is particularly true when the errors in the observed intensities are significant, 
as the first-pass solutions exactly reproduce the input data and there is no 
guarantee that a set of noisy data intensities can be exactly reproduced by a 
purely positive DEM. To mitigate this issue, we once again seek to minimize the 
total squared emission measure, but we now allow the DEM to deviate from the 
input data at a specified level. This is implemented by seeking a new DEM 
which minimizes the sum 

where £{T) is the DEM described above, and A is a regularization parameter 
chosen to enforce the desired threshold, Xo- Using our basic DEM solution, 
we replace with a set of corresponding regularization corrections to the data 
values, Ali, so that ej = A^^{Ik + /\Ik)hk- 

[S{T)]'dT = ^ ^ + A ^(/, + A/,)^(4 + Ah). (4) 

This nonlinear system of equations may be solved for the data corrections, A/|, 
and regularization parameter A, satisfying the desired threshold y^^Jn_a_ 



varie ty of ways; we have found that a standard bisection search (see iPress et al. 



2OO2I, for example) gives acceptable performance. The new regularized DEM is 



simply the first-pass inversion applied to the new corrected data values. 

Regularized solutio n of the DEM problem has recently been discussed by 



Hannah & Kontan (|2012l ), but our regularization is simpler and faster owing to 
the basis set used for the DEM. Since DEMs constructed from the instrument 
response functions have only a small number of basis elements, however, they 
remain liable to producing negative emission in cases with sharp features. The 
solution to this problem is discussed next. 
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2.3. Enforcing Non-Negativity 

We remove the remaining negative emission via the following iterative 
process: 

1. Zero the negative EM in the current DEM, 5^"^, to create a new DEM, 
£^^\ At the zeroth iteration, this is the regularized DEM from Section |Z2] 

2. Compute the data intensities, = J £^\T)Rj{T)dT, corresponding to 

p(n) 
^+ ■ 

3. Take the difference between and the original Ij, A/+ = — Ij. 

4. Compute correction DEM coefficients, Aej = Ylik^'jk^^k hk- ^Jk is 
computed using an SVD at this step, and we enforce a relatively strong 
minimum condition number (~ 0.1) to reduce ringing in the correction. 

5. Subtract the corresponding DEM corrections, AS^"-^ from £^^\ By 
construction, this restores = — AE'^'"'^ to agreement with the 
data, but reintroduces some negative emission. 

6. Repeat from steplDuntil matches li to within the desired x^- 

After we zero the negative EM in the initial DEM (i.e., at steplUof the zeroth 
iteration), the DEM is no longer represented by a basis of instrument response 
functions, but rather by a continuous function (see step |2] of the iteration). In 
practice, we choose to express it using an intermediate basis of closely spaced, 
narrow functions - usually triangle functions. 

To further speed convergence of the iteration, we attempt a linear 
extrapolation at each step of the iteration. The extrapolation steps move the 
DEM vector further along the direction of the current iteration step, and are 
only accepted when they improve x^- 

We find that the optimal regularization strength, and minimum condition 
number strength in step |4] of the iteration, vary with input DEM and data 
quality, so we try multiple values of these parameters, beginning with a light 
regularization and small minimum condition number threshold and moving 
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to stronger regularization and larger minimum condition number threshold. 
Typical values of these parameters are [0.9, 0.5, O.Olj^ and [0.01, 0.05, 0.1], for 
the regularization strengths and minimum condition numbers, respectively. In 
broad terms, the operation of the algorithm can then be described as follows: 

• Apply the light (e.g., 0.9) regularization to the data and compute the 
corresponding first-pass DEM. 

• Zero the negative emission in the first-pass DEM and compute the 
of the resulting data with respect to the initial data. If exceeds the 
desired threshold, attempt to iterate away the negative emission, using the 
smallest minimum condition threshold (e.g., 0.01), until the threshold 
is reached. 

• If the iteration takes too long to reach the desired x^ threshold, retry 
the first two steps with the next strongest regularization strengths and 
minimum conditions (e.g., 0.5 and 0.05, respectively). Repeat until the 
X^ threshold is reached, or each pair of regularization strengths and 
minimum conditions have been tried. 

Readers who would like more implementation details are encour- 
aged to examine our code, which we have made publicly available at 
|http :/"/ solar . physics . montana . edu/plowm an/f irdems .tgz , We also 
intend to submit the code to SolarSoft in the near future. 



3. Fidelity & Performance 

We have tested our DEM using a variety of example cases. The results 
are given in figures on the following pages, which show recovered DEMs, x^ 
agreement with the data, and computational time. Please note the following: 



4 



We define the regularization parameters to be p-values of the x^ distribution 



(see, for instance. Chapter 15 of iPress et al.ll2002l, where they are referred to as 



Q, and P refers to their complement). With a regularization parameter of 0.9, 
for instance, the x^ of the regularized data with respect to the original data will 
correspond to a p-value of 0.9. 
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• All reported have been reduced by dividing by the number of 
instrument channels being fit. 

• Unless otherwise noted, we use a one second exposure time for AIA, and 
30 seconds for EIS. EIS pixels were assumed to be sampled to 2 square 
arcseconds. 

• AIA errors are assumed to be from read and shot noise alone, while EIS 
errors assume an additional 10% error from other sources. 

• We use the 94, 131, 171, 193, 211, and 335 A AIA channels throughout. 

• EIS emissivities were computed assuming a density of 10^ cm^'^, unless 
otherwise noted. 

• AIA response functions were computed by calling aia_get_response ( / temp, 
/dn, /evenorm). 

• Computation times are for a single thread running on a 3.2GHz Intel Xeon 
processor, in IDL. 

Wherever possible, we use \ogiQ{T) as a temperature variable rather than 
T (in Kelvin) itself. We believe this is the more natural parameter for the 
temperature, since we are interested in a temperature range spanning multiple 
orders of magnitude (~ 10^ — 10'' Kelvin). Similarly, our DEMs are scaled per 
unit logioiT), rather than per unit T. We believe this is the natural scaling 
when the DEMs are represented as functions of \og^Q{T), because the area under 
the DEM curves is then the emission measure. Rescaling to unit T may be 
accomplished by dividing by T In (10). 

With few exceptions, we find that we are able to recover the test cases with 
good (Reduced of orie or two), that the recovered DEMs are a reasonable 
qualitative match to the input DEMs. Typical times for AIA DEMs are 
approximately one millisecond, with some cases taking under 0.1 millisecond. 
The test cases are as follows: 

• AIA inversion of Log-normal DEMs with widths of 0.2 at selected 
temperatures (Fig. [1). 
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EIS inversion of Log-normal DEMs with widths of 0.2 at selected 
temperatures (Fig. |2]). 

AIA inversion of Log-normal DEMs at temperatures ranging from 10^-^ to 
10^ Kelvin, for widths of 0.2 (Fig. 13) and 0.3 (Fig. HJ). 



EIS in version (using the 24 lines chosen by Warren et al JWarren et al. 



(|2011l) ) of Log-normal DEMs at temperatures ranging from 10^-^ to 10^ ° 



Kelvin, for widths of 0.2 (Fig. jS)) and 0.3 (Fig. ^. 

AIA inversion of DEMs produced by summing five Log-normal DEMs 
with randomly chosen centers, widths, and amplitudes (Fig. [7|. 



• EIS Active region DEM from I Warren et al.l (|201lh (Fig. H]). A density of 
]^g9.5 is assumed, to match their emissivities. 

We find that narrow (compared with the temperature response functions in 
question) DEMs are the most difficult for our method to reconstruct, in terms of 
obtaining good (i-e., Xr ~ 1) without negative emission. This can be seen by 
comparing Figures |3] and |4] (for AIA), or Figures |5] and [6] (for EIS). For AIA, there 
is somewhat more difficulty in recovering narrow DEMs at temperatures above 
10^^ Kelvin, but even in that case, acceptable were achieved for the majority 
of noise realizations. Despite their relatively poor x^, these reconstructed 
DEMs are well-behaved and localized at the center of the injected DEM; the 
median logio(T) is recovered to within 0.1 between logio(T) ^ 5.8 and 7.0. The 
difficulty in achieving good is absent in cases where there the emission is not 
concentrated near a single temperature, as can be seen in figure [71 In all cases 
considered, the average time to compute a DEM was under 10 milliseconds, 
and less than 0.1 millisecond for some cases. 

In the case of EIS, temperatures greater than 10^-^ Kelvin, above the 
highest temperature peak of the spectral lines used, are the most challenging. 
The median temperatures are accurately (i.e., 6\ogiQ{T) < 0.1) recovered 
over a range of \ogiQ(T), 5.6 ... 6.8. Due the larger number of channels, the 
computational time is somewhat longer, at ~ 10 milliseconds for a narrow 
log-normal distribution. 

Figure [Zlcompares o ur DEM inversions (top set) with those of the optimized 



Hannah & Kontaii (|2012h (bottom set) 'demmap' code. We ran the their code 
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Fig. 1. — DEM inversions of Log-normal simulated DEMs of width 0.2 using the 
six AIA EUV channels. The input DEM is shown by the solid line, while each 
dotted line is a DEM inversion with randomly chosen read and shot noise. 
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without its positivity constraint option, finding better were achieved by 
simply zeroing the negative emission. The tests DEMs are randomly generated 
multimodal distributions, as observed by AIA. Both codes were run in a single 
thread on the 3.2 GHz processor mentioned above. The results are qualitatively 
quite similar, although the Hannah & Kontar method appears to produce some 
spurious high-temperature emission. For these test cases, our DEM method 
is faster at ~ 10"'^ seconds for most DEMs compared with 3.7 x 10^^ seconds 
(without positivity constraint) for the Hannah & Kontar DEM. Our method also 
appears to produce reasonable with zero negative emission more consistently 
than the Hannah & Kontar DEM inversions. 

Figure [8] compare s our DEM results w ith the PINTof ALE MCMC results 
reported in Figure 4 of I Warren et al.l (|201l[) . These results include all 24 EIS lines 
used by Warren et al, with their factor of 1.7 tweak to the Mg intensities, along 
with the X RT Open/ Al- thick fi lter. Our DEM curve is considerably different 
than that of Warren et al. (|2nill) . In particular, we do not find a peak in the DEM 
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Fig. 2. — Same as figure [H but for the 24 EIS lines used in lWarren et al.l (|201lh 
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Fig. 3. — AIA Response to Log-normal DEMs of width 0.2 and total EM 
5.0 X lO^^cm"^ at temperatures from 5.5 to 7.0 dex. Left: the recovered DEMs 
- each vertical slice of this plot is a DEM like one of the output curves of Fig. 
[H at the temperature of the corresponding x axis value. The solid lines on left 
show emission measure weighted median temperature (EMWMT). Center: 
percentiles resulting from repeated MC trials of the read and shot noise, at each 
temperature. Right: The average time per DEM at each temperature. 
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Fig. 4. — Same as Fig. [3l but for Log-normal DEMs of width 0.3. 
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Fig. 5. — EIS Response to Log-normal DEMs of width 0.2 and total EM 5.0 x 
lO^^cm"^ at temperatures from 5.5 to 7.0 dex. Solid lines on left shows emission 
measure weighted median temperature (EMWMT). The 24 spectral lines from 
Warren et al. (2011) were used. 
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at \ogiQ{T) ^ 6.6. Despite this difference, we obtain a very reaso nable of 



1.16, w hich is fully consistent with the data. This suggests that the IWarren et al, 



(|201l[) reconstruction may not represent the full range of DEMs consistent with 
their data. 



3.1. Example Solar DEM Analysis 

Finally, we show example DEM analysis of solar data. Figure |9] shows 
emission weighted median temperature (EMWMT) and total emission measure 
for an active region observed by EIS and AIA. The data cover the area of the 
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Fig. 6. — EIS Response to Log-normal DEMs of width 0.3 and total EM 5.0 x 
lO^^cm"^ at temperatures from 5.5 to 7.0 dex. Solid lines on left shows emission 
measure weighted median temperature (EMWMT). The 24 spectral lines from 
Warren et al.l (|201l|) were used. 
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large raster taken by EIS at around 01:30 on April 19, 2011. They are centered 
~ 350 arcseconds above disk center, with an approximately 230 x 510 arcsecond 
field of view. The EIS data used fits to a set of five iron lines: Fe IX 188.497, Fe X 
184.537, Fe XII 195.119, Fe XV 284.163, and Fe XVI 262.976. 

The average time per AIA DEM was approximately 0.13 milliseconds, and 
we achieve reasonable for over 95% of the AIA DEMs (the 95th percentile 

for EIS is 4.7). The active region plots in Figure |9] are qualitatively quite 
similar, suggesting that the results give useful insight into the underlying solar 
temperature distribution. 

We also show DEMs along the loop segment indicated by the dashed line 
in Figure m Figure [10] shows DEMs along the length of the loop segment (solid 
lines in Figure [TTj) , as well as equivalent sets of DEMs offset to either side of the 
loop (dotted lines in Figure [TTj) . Once again, the DEMs found by AIA and by 
EIS are reasonably consistent, and both EIS and AIA inversions achieved good 
for over 95% of the points along the loop. 

The DEMs shown along the loop axis are not made from background 
subtracted data, and they show little difference from the offset DEMs, 
which indicates that the loop emission has been swamped out be the bright 
background. In an upcoming work, we compute background subtracted DEMs 
for this loop and compare its temperature and density profile to a set of analytic 
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strand heating models. 



4, Conclusions 

We have demonstrated a method for fast reconstruction of DEM 
distributions using coronal data from instruments such as EIS and AlA. This 
DEM method achieves reduced chi-squared of order unity with no negative 
emission in all but a few test cases. The most difficult test cases are narrow 
DEMs at high (> 10^'^ Kelvin) temperatures for AIA, and temperature regions 
with little spectral coverage for EIS. Even for the high temperature AIA 
cases, we achieve reasonable for the majority (~ 75%) of noise realizations. 
Qualitatively, the reconstructed DEMs match the input DEMs well, although 
the ability to recover finer details of the input DEMs is inherently limited. 
The data, particularly for AIA, do not constrain the fine details of the DEMs. 
When interpreting DEMs, great care must be exercised to determine whether 
or not the features of interest are genuine aspects of the coronal temperature 
distribution. 

Our DEM method is fast enough to study the dynamics of the coronal 
temperature distributions in real time. For AIA data, the worst cases 
considered take under 10 milliseconds, and some test cases execute in less 
than 0.1 milliseconds. For the whole disk, at full AIA resolution using all 
six coronal EUV channels, we expect less than ~ 1 hour with the code as 
currently implemented. We plan to convert the computationally intensive 
parts of the code to C and rewrite them to take advantage of multithreading, 
which should offer a factor of ~ 100 performance gain on an eight core 
workstation. This would reduce the time for a set of full-resolution, full-disk 
AIA DEMs to under ~ 1 minute, and achieve the objective of matching 
the AIA observing rate in real time. The software is available online at 
Ihttp : //solar .physics .montana . edu/pl owman/f irdems . tgz l and 
will be submitted to SolarSoft in the near future. 

We also applied our DEM reconstruction to solar data, an active region 
observed by EIS and AIA on April 19, 2011. We found no difficulty achieving 
reduced of order unity with no negative emission, and the average time per 
DEM was approximately 0.13 milliseconds, or 44 seconds for the entire active 
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region at full AIA resolution. We find a relatively hot active region core with 
median temperature of around three million Kelvin, surrounded by cooler 
emission with median temperature of around two million Kelvin. We also 
plotted DEMs as a function of length along a coronal loop segment visible in the 
active region. These DEMs are dominated by a relatively bright background, 
making it unclear what emission is associated with the loop; background 
subtraction of the data will be necessary to isolate the loop emission from its 
surroundings. 



A. Equality of SVD-derived and Instrument Response Basis DEMs 

The fundamental equation we want to invert is the following: 

I. = J Rj{T)£{T)dT (Al) 

Without loss of generality, we can represent this integral in matrix form: 

/ = R^;, (A2) 

where the elements of the n x m (for ric channels and Ut temperature bins) 
response matrix R are Rj{Tk), and the elements of the rif-element DEM vector 
E are £{Tj)AT. The singular value decomposition of R is 

R = USV^, (A3) 

Where U is an ric x Uc unitary matrix, V is an rit x rit unitary matrix, and S is 
the ric X ut diagonal matrix of the singular values of R. The SVD can be used to 
estimate E by computing the pseudo-inverse of of R, R+ = VS^U^: 

^SVD _ (^4) 

We can also choose to represent the DEM using a basis of instrument response 
functions: 

E^^ = R^e (A5) 

The ric-element vector of coefficients e is obtained from I by plugging E^^ into 
equation I A2I and inverting: 

e = [RR^]-ij (A6) 
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Expanding [RR^]~^ in terms of the SVD of R, we find that 

[RR^]-^ = [USV^VS^U^]-^ = [USS^U^]-^ = U(S+)^S+U^. (A7) 
We then obtain 

E^^ = RT[RR^]-^/ = VS^U^(U(S+)^S+U^)J = VS+U^J = R+J. (A8) 
Therefore, E^^ = E^^^, as claimed. 
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Fig. 7. — AIA inversion of DEMs produced by summing five Log-normal DEMs 
with randomly chosen centers, widths, and amplitudes. Center locations are 
uniformly distributed between 10^'^^ and 10^, widths between 0.1 and 0.3 and 
total EM between 5 x lO^^cm"^ and 5 x 10^^cm~^. Top : Results from our fast DEM 



Method. Bottom: Same for iHannah & Kontarl (|2012l) regularized DEM method, 
without positivity constraint. 



0.080582 ms per DEM; 95% Chi St]uarccl= 1 .30 




Log,„(T) (Kelvini 
0.084198 ms per DEM; 95% Chi Squarecl= 1 .09 




Log|„(T) (Kelvin) 
_^0.1956i m:i per DEM; 95% Clii Squai-ed= 2.09 





Log,o(T) (Kelvin) 
,s per DEM; 95% Chi Squiirec]= 2.08 




Logi„(T) IKelvin) 
0.079175 ms per DEM; 95% Chi Squiircd= 1.17 




Log|„(T) IKelvin) 

0.083583 ms per DEM; 95% Chi Squarcf]= 1 .04 
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Fig. 8. — Comparison of our DEM results (top) with lWarren et alj (|201ll) (Figure 
4) MCMC DEM results (bottom). We include their factor of 1.7 tweak to the 
coronal Mg abundance. The DEMs matched the data with of order unity. 
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AIA Median Temperatures EIS Median Temperatures Color Map 




Fig. 9.— Left: EMWMT (hue) and Total EM (intensity) from AIA DEM inversion 
of Active region (covered by EIS fov on April 19, 2011). Middle: Same for EIS. 
Right: color scale for middle and left plots. 



AIA OEMs (Offset Left from Loop) AIA Loop DEMs AIA DEMs (Offset Right rniiii Loo]i] 




20 40 60 80 20 40 60 80 20 40 60 80 

Lengiti (arcseeonds) Lengili (arcseconds) Lengih (arcsefouds) 



Fig. 10. — Top: Loop DEM from AIA inversion of Active region (Center), offset 5 
arcseconds to the left (left), 5 arcseconds to the right (right). Loop area indicated 
by dashed line in Figure HI or in more detail in Figure [TTl Bottom: Same for EIS. 
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Fig. 11. — Detailed, straightened AIA images of loop outlined in Figure HI for 
comparison with Figure [lOl Loop trace indicated by solid line, while ±5 arcsec- 
ond offset background traces are indicated by dotted lines. 
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